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1.  Introduction 


This  technical  note  explores  the  design  of  the  monolithic  neodymium  (Nd):  yttrium  aluminum 
garnet  (YAG)  laser  used  in  the  diode-pumped  laser  ignition  system  (DPLIS).  Emphasis  is  placed 
on  the  divergence  of  the  output  beam  and  the  dependence  of  the  output  power  on  misalignment 
of  the  mirrors  polished  on  the  two  ends  of  the  rod.  The  design  parameters  include  the  radius  of 
curvature  of  each  mirror  and  the  length  of  the  rod.  The  DPLIS  is  a  candidate  to  replace  the 
primer-based  ignition  of  propellant  in  a  155-mm  howitzer.  The  DPLIS  currently  meets  all  the 
required  specifications,  but  additional  fluence  would  be  desirable  under  conditions  where  the 
distance  to  the  target  is  as  large  as  30  in  (0.762  m)  (1).  Hence,  the  importance  of  minimizing  the 
divergence. 


2.  Fundamental  Mode  Properties 


The  rod  geometry  and  its  relationship  to  the  target  are  shown  in  figure  1.  R  i  and  Ri_  are  the  radii 
of  curvature  of  the  two  mirrors.  L  is  the  length  of  the  rod  and  T  is  the  distance  to  the  target.  Not 
counting  the  cladding  surrounding  the  Nd-doped  portion,  the  diameter  of  the  rod  is  2 r.  A  high 
reflectivity  coating  is  applied  to  Mi  and  a  partially  reflecting  coating  applied  to  M2,  the  output 
coupler. 


Figure  1.  Geometry  of  the  laser  rod,  cladding,  and  target. 

We  apply  the  usual  analysis  to  determine  the  fundamental  (TEM0o)  mode  inside  and  outside  the 
cavity  (5).  The  fundamental  mode  is  defined  if  one  knows  the  complex  beam  radius,  q,  at  any 
point  along  the  axis  of  propagation. 


1  _  1  a 

q  p  mux)2 


1 


where 


A  is  the  wavelength.  For  the  rest  of  this  technical  note,  A  —  1064  nm. 
p  is  the  radius  of  curvature  (r.o.c.)  of  the  wavefront 

n  is  the  index  of  refraction  of  the  medium.  For  YAG  at  1064  nm,  n  =  1.82. 

( o  is  the  spot  size.  At  a  distance  w  from  the  axis,  the  intensity  has  decreased  by  1/e  . 

The  complex  beam  radius  incident  upon  Mi  is  given  by  the  solution  of 

Aq 2  T  B  1 

q2  —  — - —  or  — 

Cq2  +  D  q2 

The  latter  form  is  convenient,  given  the  relation  of  q  to  p  and  co. 
rod,  the  roundtrip  matrix  is  given  by 

(A  Bill  LU  1  0W1  LU  1  0\ 

VC  D>  VO  V\-2/R,  ljlo  1/  \—2/R2  1) 

When  R2  =  oo,  a  waist  forms  at  Mi  with  a  spot  size  ah  given  by 


C  +  D  /  CJ2 

= - (2 

A  +  B  /  q2 

When  there’s  no  lensing  in  the 


0>2  = 


-  ^0  -  ^  n  v  ^(^1 


The  radius  of  curvature  of  the  beam  matches  that  of  the  mirror  at  each  end  of  the  cavity.  As  an 
example,  we  calculate  the  fundamental  mode  for  the  case  =  5  m,  L  —  10  cm,  R2  =  00.  The 
first  column  of  figure  2  shows  the  spot  size  and  radius  of  curvature  inside  the  cavity.  In  the 
second  column,  the  same  quantities  are  shown  between  the  laser  and  the  target.  The  dashed  lines 
indicate  the  asymptotic  dependence.  The  Matlab©  program  that  generated  the  data  is  reproduced 
in  the  appendix. 
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Figure  2.  TEM0o  mode  of  the  resonator,  when  K]  =  5m ,  L  —  10  cm,  R2  —  °°.  Spot  size 

(top  row)  and  radius  of  curvature  (bottom)  inside  the  laser  (left  column)  and  outside 
the  laser  (right).  The  target  is  located  T  =  0.762  m  from  M2. 

Outside  the  cavity,  in  air,  the  divergence  of  the  fundamental  mode  is  given  by 

A  An 

do  = - =  r  ,  (5 

Table  1  shows  co i,  <xh,  and  6b  for  various  choices  of  L  and  R\,  for  the  case  R2  —  00 .  R\  >  10L  in 
all  cases.  Two  points  can  be  made:  (1)  The  spot  sizes  at  the  two  mirrors  are  very  close  to  one 
another,  because  the  rod  is  short  compared  to  R\  and  R2  in  all  cases  and  (2)  as  R\  increases,  0(} 
decreases. 


3 


Table  1.  Spot  sizes,  divergences,  and  alignment  sensitivity  for  specific  values  of  L,  Rh  and  R2. 


L 

(m) 

Ru 

HR 

(m) 

Ri, 

OC 

(m) 

®1 

(mm) 

co2 

(mm) 

(mrad) 
half  Z 

M 

0m 

(mrad) 
half  Z 

emr 

(mm) 

<bt 

(mm) 

D 

(xlO4) 

(co  tD) 
(xlO3) 

0.1 

1 

00 

0.2491 

0.2363 

1.43 

10.6 

15.17 

11.56 

11.79 

0.596 

70.3 

ff 

2 

ff 

0.2922 

0.2848 

1.19 

8.78 

10.44 

7.95 

8.29 

1.02 

84.2 

ff 

5 

ff 

0.3646 

0.3609 

0.938 

6.93 

6.50 

4.95 

5.48 

2.03 

111 

ff 

10 

ff 

0.4325 

0.4303 

0.787 

5.81 

4.57 

3.48 

4.20 

3.43 

144 

ff 

20 

ff 

0.5136 

0.5124 

0.661 

4.88 

3.23 

2.46 

3.39 

5.77 

196 

ff 

40 

ff 

0.6104 

0.6097 

0.556 

4.10 

2.28 

1.74 

2.91 

9.71 

283 

.075 

1 

ff 

0.2302 

0.2214 

1.53 

11.3 

17.27 

13.16 

13.37 

0.645 

86.2 

ff 

2 

ff 

0.2710 

0.2649 

1.27 

9.40 

11.97 

9.12 

9.42 

1.09 

103 

ff 

5 

ff 

0.3389 

0.3363 

1.01 

7.43 

7.49 

5.70 

6.17 

2.19 

135 

ff 

10 

ff 

0.4022 

0.4007 

0.845 

6.24 

5.27 

4.02 

4.65 

3.69 

171 

ff 

20 

ff 

0.4778 

0.4770 

0.710 

5.24 

3.72 

2.84 

3.68 

6.20 

228 

ff 

40 

ff 

0.5680 

0.5675 

0.597 

4.41 

2.63 

2.00 

3.08 

10.4 

322 

3.  Multimode  Properties 


For  typical  values  of  R i  and  Ri,  a  5-mm-diameter  rod  (r  =  2.5  mm)  will  lase  with  multiple 
transverse  modes.  The  laser  can  be  forced  to  lase  in  the  fundamental  mode,  however,  by 
introducing  an  aperture  -20%  larger  than  the  TEMoo  spot  size.  Theoretically,  one  can  also 
increase  an  and  ah  to  the  point  where  the  finite  rod  diameter  suppresses  the  higher-order  modes. 
However,  this  is  impractical  because  it  would  require  R1  S  5  km. 

The  far-field  divergence  of  a  higher-order  mode  is  given  by 

em  =  Vm  +  1/2  0O  (6) 

where  m  is  the  transversal  order  of  the  mode.  According  to  reference  4,  if  an  aperture  is  placed 
inside  the  resonator,  and  the  rod  is  uniformly  excited,  the  number  of  oscillating  modes  can  be 
estimated  by 


m  +  1/2  =  (a/u0a)2  (7) 

where 

a  is  the  radius  of  the  aperture.  For  the  remainder  of  this  technical  note,  a  =  r  —  2.5  mm. 

a>0a  is  the  spot  size  of  the  fundamental  mode  at  the  position  of  the  aperture. 

In  our  case,  rather  than  a  conventional  aperture  at  one  fixed  position,  the  Nd:YAG  rod 
constitutes  an  aperture  along  the  whole  length  of  the  rod,  because  the  gain  is  zero  everywhere 
outside  the  rod  radius.  However,  the  aperture  is  “soft”  because  the  mirror  diameter  is  equal  to 
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that  of  the  rod  plus  the  cladding,  so  light  travelling  outside  the  rod  can  remain  in  the  cavity, 
albeit  with  lower  gain. 

When  R2  —  oo,  the  multimode  divergence  is  given  by 

nr  Aa 

dm  —  i  „  —  2  " 

Vz.(/?1  -L) 

The  divergence  has  a  minimum  of  dm  =  1.82  mrad  for 
half  of  a  confocal  resonator. 

Note  that  if  m  modes  are  superimposed,  with  random  phases,  the  beam  would  have  a  beam 
quality  of  M2  —  m  +  1/2  in  each  transverse  direction.  The  divergence  and  the  spot  size  are 
each  M  times  larger  than  the  values  for  the  fundamental  mode;  therefore,  the  net  beam  quality  is 
M 2. 

It  is  not  clear  why  the  divergence  of  the  entire  beam  should  be  equal  to  that  of  the  highest-order 
mode.  Intuition  would  have  it  equal  to  some  kind  of  average  divergence  of  all  the  modes. 

As  a  multimode  beam  propagates  a  distance  T,  the  spot  size  increases  approximately  as  (5) 

6>t  =  wm2  +  (dmT)2  (9) 

where 

coT  is  the  radius  of  the  spot  that  contains  87.5%  of  the  power  at  the  target,  and 

com2  is  the  multimode  spot  size  at  the  beam  waist,  which  is  located  at  M2  when  R2  =  oo.  Because 
the  multimode  spot  is  not  Gaussian,  it  is  not  obvious  how  to  define  the  spot  size,  so  we  use  the 
87.5%  criterion.  The  spot  size  of  a  TEMoo  beam  contains  87.5%  of  the  power.  For  a  uniform 
intensity  profile  inside  the  2.5-mm-radius  rod,  a  spot  size  of  com2  =  2.34  mm  would  contain 
87.5%  of  the  power. 

0m  is  the  multimode  divergence. 

T  is  the  distance  to  the  target. 

Table  1  shows  the  calculated  values  for  M,  6m,  0mT,  and  oT  for  various  choices  of  L  and  R\,  for 
the  case  a  =  2.5  mm  and  T  =  0.762  m.  Four  points  can  be  made: 

1.  As  R2  increases,  both  M  and  0O  decrease,  therefore  0m  decreases. 

2.  For  the  shorter  values  of  R\,  a> r  is  dominated  by  the  divergence,  i.e.,  coT  ~  9mT. 

3.  For  the  larger  values  of  Ri,  e.g.,  Rt  =  40  m  and  L  =  0.1  m,  the  divergence  contributes  only 
about  1/3  to  the  spot  size;  the  rest  is  due  to  <wm2. 

4.  The  shorter  rod  length  results  in  a  larger  car. 


a 

do  — 

O)0 


(8) 


=  2 L.  This  geometry  is  equivalent  to 
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2 

The  6.5  mrad  (half  angle  l/e“)  in  row  three  agrees  well  with  the  9-mrad  divergence  (full  angle 
50%  power)  measured  experimentally*.  Were  the  beam  Gaussian,  the  latter  would  be  equivalent 
to  a  7.64-mrad  half  angle  1/e  divergence. 


The  contour  plots  in  figure  3  show  co2  (for  the  fundamental  mode),  M,  dm,  and  cox  as  a  function  of 
R\  and  L,  for  R2  —  o°.  The  spot  size  at  the  target  clearly  decreases  for  larger  values  of  R\  and  L. 
However,  we  expect  the  alignment  sensitivity  to  increase  as  well. 


to2  spot  size  at  M2  (mm) 


M 


E, 

E 


L  (no) 

coT  spot  size  at  target  (mm) 


E 


L  (m) 


Figure  3.  Contour  plots  of  various  quantities  as  a  function  of  R\  and  L,  for  R2  — 

(upper  left)  ah,  spot  size  of  the  fundamental  mode  at  M2;  (upper  right)  M,  i.e.,  the 
ratio  of  rod  radius  to  spot  size  at  M2;  (lower  left)  9m,  divergence  of  multimode  beam 
(half  angle  of  cone  that  contains  87.5%  of  power);  and  (lower  right)  <rur,  spot  size  at 
target  0.762  m  away  from  M2. 

Next,  we  consider  the  impact  of  a  finite  positive  radius  of  curvature  on  Mi.  For  the  case 
L  —  0.1  m,  the  spot  sizes  and  divergence  are  shown  in  figure  4. 


* 

The  output  coupler  had  a  reflectivity  of  50%  (6). 
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gi2  spot  size  at  M2  (mm) 


em  (mrad) 


toT  spot  size  at  target  (mm) 


Figure  4.  Same  quantities  as  in  figure  3,  plotted  as  a  function  of  R\  and  R2,  for  L  —  0.1  m. 


4.  Alignment  Sensitivity 


The  alignment  sensitivity  for  a  fundamental  mode  two-mirror  cavity  with  apertures  at  each 
mirror  has  been  analyzed,  and  the  analysis  verified  with  experiment  (7).  Assuming  small 
misalignments,  i.e.,  small  losses,  and  assuming  that  the  apertures  at  each  end  of  the  rod  are  1.2x 
the  spot  size,  the  misalignment  sensitivity  of  Mi  is  given  by 


1  +  9i9i 

(1  -9i92)3/2 


(10) 


where  gt  —  1  —  L/Ri.  If  M;  is  tilted  by  an  angle  of  1/D;,  additional  losses  of  -10%  are  incurred 
(4).  The  combined  sensitivity  of  the  two-mirror  cavity  is  defined  as 


D2  =  D{  +  D2 


(11) 
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The  alignment  sensitivity  would  be  at  minimum  for  a  confocal  geometry  (not  shown  in  table  1). 
A  disadvantage  of  the  symmetric  confocal  geometry  is  the  marginal  stability  in  the  ABCD  matrix 
sense.  Very  small  changes  to  Ru  /A,  or  L  can  make  the  cavity  unstable. 

Equations  10  and  1 1  apply  to  a  fundamental  mode  cavity.  We  make  the  assumption  that  they 
also  apply  to  a  multimode  cavity  if  the  aperture  is  equal  to  1.2x  the  multimode  spot  size,  which  is 
approximately  the  case.  The  resulting  values  of  D  are  shown  in  table  1.  The  trend  should  be 
meaningful  even  if  the  values  are  not  quantitatively  correct. 


The  top  row  of  figure  5  are  plots  of  D\  and  D2  as  a  function  of  R\  and  L,  for  R2  —  00.  It  is  clear 
that  to  reduce  the  sensitivity,  one  would  want  to  decrease  R\.  In  the  regime  we  are  considering, 

L  has  less  of  an  impact  on  D.  The  contour  D  —  20,000  corresponds  to  a  misalignment  of 
50  prad  introducing  a  loss  of  10%.  Even  though  R\  and  R2  are  quite  different,  g  1  and  go  are  both 
close  to  one;  therefore,  D\  is  nearly  the  same  as  Do.  Figure  5c  shows  the  misalignment  angle  that 
would  introduce  an  additional  loss  of  10%,  based  on  the  combined  sensitivity. 


sensitivity  D1 


sensitivity  D, 


10 

8 

g,  6 
E  4 
2 


2000Q 

15000 

-10000 


0.05  0.1  0.15 

L  (m) 

10%  loss  angle  1/D  prad 


cr 


cr 


L(m) 

cot*D 


L  (m) 


Figure  5.  Contour  plots  of  various  quantities  as  a  function  of  Rl  and  L,  for  R2  = 
(top  left)  Mirror  misalignment  sensitivity  D\,  (top  right)  D2,  (lower  left) 
1  ID,  and  (lower  right)  o>rD . 
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We  take  the  product  coTD  to  be  a  figure  of  demerit,  i.e.,  the  larger  the  product  the  less  desirable 
the  cavity.  The  values  of  0)\D  for  specific  values  of  R\  are  shown  in  table  1.  The  product  is 
shown  in  figure  5  (lower  left)  as  a  function  of  R\  and  L.  The  contour  labeled  100  corresponds, 
for  example,  to  a  divergence  of  1  mrad  and  an  alignment  tolerance  of  10  prad. 


The  same  quantities  are  plotted  in  figure  6,  as  a  function  of  R\  and  R2,  for  the  case  L  —  0.1  m. 


E 

EE 


E 

EE 


sensitivity  D1 


10%  loss  angle  1/D  grad 


2  4  6  8  10 


R2  (m) 
toT  *  D 


Figure  6.  Contour  plots  of  the  same  quantities  as  in  figure  5,  plotted  as  a  function  of 
and  R2,  for  L  —  0.1  m. 


5.  Discussion 


The  three  design  parameters  that  can  easily  be  adjusted  are  L,  Ru  and  R2.  According  to  the 
above  results,  a  longer  L  is  desirable;  however,  L  is  limited  to  10  cm  by  the  space  available  in  the 
breech.  It  is  also  clear  from  the  preceding  analysis  that  larger  mirror  radii  of  curvature  yield  a 
smaller  spot  on  the  target,  and  thus  higher  fluence  for  a  given  pulse  energy.  However,  a  choice 
of  optimum  R\  and  R2  has  to  be  informed  by  the  difficulty  (or  cost)  in  meeting  the  required 
parallelism  tolerance,  which  increases  with  R.  Another  avenue  to  pursue  is  an  optimization  of 
the  output  coupler  reflectivity.  Adding  a  lens  next  to  the  pressure  window  might  also  be  helpful. 
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Making  the  window  into  a  positive  lens  by  curving  one  surface  is  also  being  explored  (i). 


6.  Conclusion 


The  analysis  presented  in  this  report  quantifies  the  known  tradeoffs  associated  with  varying  the 
radius  of  curvature  of  the  mirrors  polished  on  the  rod  in  the  DPLIS.  Increases  in  radius  of 
curvature  decrease  the  divergence  of  the  output  beam,  which  is  good,  but  they  also  make  more 
stringent  the  parallelism  of  the  two  mirrors,  which  adds  to  the  cost  of  fabrication. 
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Appendix.  Matlab  Program  that  Generated  the  Data 


The  following  is  the  Matlab  program  that  generated  the  data. 

function  [ zv, wv, rhov, w2 , thetaO , M, D1 , D2 , spotT]  =  wrho (R1 , L, R2 , T, printf ig) 

%  outputs  z, w, rho, thetaO,  M,  D,  spot  size  at  target  given  R1,L,R2,T 
%  when  printfig  is  set  to  1,  a  figure  is  displayed  showing  the  results 
global  wl  n  r 

fprintf('Rl=  %6.2e  m\n',Rl);  % 

fprintf('L  =  %6.2e  m\n',L);  % 

fprintf('R2=  %6.2e  m\n',R2);  % 

fprintf('T  =  %6.2e  m\n',T);  % 

M2  =  [1  L;  0  1]  *  [1  0,-2/Rl  1]  *  [1  L;  0  1]  *  [10;  -2/R2  1];  %  at  M2 

A2  =  M2  ( 1 , 1 )  ;  B2  =  M2  ( 1 , 2  )  ;  C2  =  M2  ( 2 , 1 )  ;  D2  =  M2  ( 2 , 2  )  ; 

fprintf('A=  %6.2e  B=  %6.2e  C=  %6.2e  D=%6 . 2e\n\n 1 , A2 , B2 , C2 , D2 ) ; 

ooq2i  =  ( (D2-A2) -sqrt ( (D2-A2) A2+4*B2*C2) ) / (2  *  B2); 

fprintf ( ' l/q2i  =  %6.2e  %6.2e  mm\n ' , real (ooq2i)  ,  imag (ooq2i)  );  % 

%  ooq2i  =  (D2-A2)/(2  *  B2)  -  li  *  sqrt(l-  (  (D2+A2 ) /2 ) A2 )  /  B2;  %  from  Kortz 

%  fprintf (' l/q2i  =  %6.2e  %6.2e  mm\n real (ooq2i) , imag (ooq2i)  );  % 

q2i  =  l/ooq2i; 

fprintf ('q2i  =  %6.2e  %6.2e  mm\n real (q2i) , imag (q2i)  );  % 

w2 i  =  w_qn ( q2 i , n ) ; 
w2  =  w2 i ; 

q2e  =  q2i  /  (  (l-n)/R2*q2i  +  n  );  %propagate  through  OC  into  air 

rho2e  =  rho_q(q2e); 

w2e  =  w_qn(q2e,l);  %  index  of  1 

thetaO  =  atan (wl  /  (pi*w2e)  ); 

M  =  r  /  w2e;  %  how  many  time  diff  ltd 

thetam  =  thetaO  *  M; 
gl  =  1-L/R1; 
g2  =  1-L/R2; 

prefactor  =  (pi*L/wl) * (l+gl*g2) / (l-gl*g2) A1 . 5; 

Dls  =  prefactor*sqrt (g2/gl)  ; 

D2s  =  prefactor*sqrt (gl/g2)  ; 

D1  =  sqrt (Dls) ; 

D2  =  sqrt (D2s) ; 

D  =  sqrt (Dls+D2s) ; 

Da  =  D; 

Ds  =  pref actor*abs (gl+g2) /sqrt (gl*g2)  ; 

Db  =  sqrt (Ds) ; 

spotTs  =  0.875*rA2  +  (thetam  *  T)A2; 
spotT  =  sqrt ( spotTs )  ; 

ziv  =  linspace (-L,  0, 100) ; 
qiv  =  q2i  +  ziv; 
wiv  =  w  qn(qiv,n); 
rhoiv  =  rho_q(qiv); 
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zev  =  linspace (0,  T,  100) ; 

qev  =  q2e  +  zev; 

wev  =  w  qn(qev,l); 

rhoev  =  rho_q(qev); 

wav  =  zev  *  wl  /  pi  /  w2e;  %  waist  asymptotic 

zv  =  cat (2 , ziv, zev) ; 

wv  =  cat (2 , wiv, wev) ; 

rhov  =  cat (2 , rhoiv, rhoev) ; 


fprintf ( ' spot  at  Ml 
fprintf('roc  at  Ml 
fprintf ( ' spot  at  M2 
fprintf ('roc  at  M2 
fprintf ( ' thetaO 
fprintf ( ' M=  %6.2e\n 
fprintf (' thetam  = 
fprintf (' thetam  T= 
fprintf ( ' spot  at  T= 
fprintf ( ' D1  =  % 

fprintf ( ' D2  =  % 

fprintf ( ' Da  =  % 

%  fprintf ('Db 
%  fprintf (' thetam  D 
fprintf ( ' spotT  *  D 
%  fprintf (' 1/ (spotT 


=  %10.3e  mm\n ', wiv ( 1 ) *1000 ) ; 

=  %6.2e  m\n ', rhoiv ( 1 )) ; 

=  %10.3e  mm\n ' , w2i*1000 ) ; 

=  %6.2e  m\n',rho2e);  % 

=  %6.3f  mrad\n ' , theta0*1000) ; 

\M);  % 

%6 . 2f \n ' , thetam* 1000 ) ; 

%6.2f  mm\n 1 , thetam*T*1000 )  ; 
%6.2f  mm\n\n ' , spotT*1000) ; 

6 . 2e  \n ' , Dl) ; 

6 . 2e  \n ' , D2 ) ; 

6 . 2e  \n ' , Da) ; 

%  6 . 2  e  \ n ' , Db ) ; 

=  %6.2e  \n '  ,  thetam*Da) ; 

=  %6.2e  \n 1 , spotT*Da) ; 

D)=  %6.2e  \n\n ' , 1 / ( spotT*Da) ) 


O, 

o 


o, 

o 


o, 

o 


o, 

o 


if  printfig==l; 

figure8xl0  (1) ; 

subplot  (3,2,1,  'Visible',  'off') ; 
annotation ( ' textbox ' , . . . 

'Position ' , get (gca, 'Position ' ) , 'String ' , . . . 

{ ' wrho . m ' ; . . . 

['Rl=  '  num2str(Rl)  'm  L=  '  num2str(L)] 

['R2=  '  num2str(R2)  '  T=  '  num2str(T)] 

['spotsize  @HR=  '  num2str(wiv(  1 ) *1000 , ' %5 . 3f ' )  'mm' 
['  @OC=  '  num2str(wev(  1 ) *1000 , ' %5 . 3f ' ) ] 

['  @target=  '  num2str (wev (end) *1000 ,' %5 . 3f ') ] 


II 

\ — 1 
Q 

'  num2str (D1 , 

'  % 6 . 2e ' ) 

]  ; 

II 

CM 

Q 

'  num2str (D2 , 

'  % 6 . 2e ' ) 

]  ; 

[  '  D  = 

'  num2str (D  , 

'  % 6 . 2e ' ) 

] ; 

}) ; 

subplot (3,2,2, 'Visible', 'off'); 
annotation ( ' textbox ' ,  .  .  . 

' Position ', get (gca,  'Position ' ) ,  .  .  . 

'String ' , { dates tr (now) ; . . . 

[ ' n=  '  num2str (n) ] ;  .  .  . 

['\lambda=  '  num2str (wl*le6)  '  um' ] ; . . . 

[ ' r=  '  num2str(r)  '  mm'];... 

[ ' \theta0=  '  num2str (theta0*1000 , ' %5 . 3f ' )  ' 

mrad '];... 

[ ' M=  '  num2str (M) ] ;  .  .  . 

[ ' \thetam  =  '  num2str (thetam* 1000 ,' %6 . 2f ')] ; 
[ ' \thetam  D  '  num2str (thetam  *  D, ' %6 . 2e ' ) ] } ) 
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subplot (3,2,3) ; 
plot ( ziv*100 , wiv*1000 ) ; 
xlabel ( 'distance  (cm) ' ) ; 
ylabel ( ' spot  size  (mm) ' ) ; 

%  v  =  axis; 

%  axis ( [  v(l)  v(2)  0.36  0.365]) 
subplot (3,2,4)  ; 

plot ( zev*100 , wev*1000 ,  -');  %, zev* 100 , wav* 1000 ,'--') ; 
xlabel ( 'distance  (cm) ' ) ; 
ylabel ( ' spot  size  (mm) ' ) ; 

subplot (3,2,5) ; 

plot (ziv*100,max (-4*Rl,min (4*R1, rhoiv) )  ) ; 

xlabel ( 'distance  (cm) ' ) ; 
ylabel ( ' r . o . c .  (m) ' ) ; 

%  v  =  axis; 

%  axis ( [  v(l)  v(2)  v(3)  0]) 

subplot (3,2,6); 

plot (zev* 100, max (-5*T,min (5*T, rhoev) ),'-');  %, zev* 100, zev, ' 
xlabel ( 'distance  (cm) ' ) ; 
ylabel ( ' r . o . c .  (m) ' ) ; 

%  v  =  axis; 

%  axis ( [  v(l)  v(2)  0  R2 ] ) ; 

end 


function  w  =  w  qn(qv,n) 

%  calculate  spot  size  given  q  vector  and  index  of  refraction 
global  wl 

w  =  sqrt (  -wl  ./  (pi  *  n  *  imag(l./qv)  )  ); 

function  rho  =  rho  q(qv) 

%  calculate  radius  of  curvature  given  q  vector 
rho  =  1./  (real ( 1 . /qv) ) ; 


global  wl  n  r 
wl  =  1.064e-6; 
n  =  1.82; 
r  =  2.5e-3; 


m 

Nd : YAG  at  1.06 
rod  radius  (m) 
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